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We describe the dynamics of gluons and quarks in a rel¬ 
ativistic nuclear collision, within the framework of classical 
mean-field transport theory, by the coupled equations for the 
classical Yang-Mills field and a collection of colored point par¬ 
ticles. The particles represent the valence quarks in the col¬ 
liding nuclei. We explore the time evolution of the gauge 
field in a numerical simulation of the collision of two Lorentz- 
boosted “nuclei” on a long three-dimensional gauge lattice. 
We report first results on soft gluon scattering and coherent 
gluon radiation. 


Experiments with relativistic heavy ions at center-of- 
mass energies reaching 100 GeV/u will soon search for 
a new phase of nuclear matter, the quark-gluon plasma 
[Q,|j. One of the theoretical challenges in this context is 
the development of a description on the basis of quan¬ 
tum chromodynamics (QCD) of the processes that may 
lead to the formation of locally equilibrated superdense 
matter in these nuclear reactions. 

The possibility of a description of inelastic gluon pro¬ 
cesses by means of the nonlinear interactions of classical 
color fields has been proposed some time ago |^,Q. Re¬ 
cently, this scenario was examined in studies of the col¬ 
lision of two transverse polarized Yang-Mills field wave 
packets on a one-dimensional gauge lattice ||]. These 
calculations showed that the interaction between local¬ 
ized classical gauge fields can lead to the excitation of 
long wavelength modes in a way that is reminiscent of 
the formation of a dense gluon plasma. 

Here we address the question to which extent soft 
gluon scattering in ultra-relativistic nuclear collisions 
contributes to the excitation of gluonic modes with non¬ 
vanishing transverse momentum and, therefore, to the 
conversion of longitudinal kinetic energy into transverse 
energy. The colliding gluon fields are generated by clas¬ 
sical color sources, representing the fast moving nuclear 
valence quarks. The inclusion of particles also permits a 
comparison with perturbative calculations of gluon radi¬ 
ation from colliding quarks . 

A set of equations describing the evolution of the phase 
space distribution of quarks and gluons in the presence 
of a mean color field, but in the absence of collisions, was 
proposed by Heinz [^|j^] . This non-Abelian generalization 
of the Vlasov equation can be considered as the contin¬ 
uum version of the dynamics of an ensemble of classical 
point particles endowed with color charge and interacting 
with a mean color field. Denoting the space-time posi¬ 
tions, momenta, and color charges of the particles by xf, 


pf and qf, respectively, where f = 1,..., A^ is the particle 
index, the equations for this dynamical system read: 
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Here g is the gauge coupling constant, are the 

structure constants of the gauge group (here taken as 
SU(2)), and F^^'' denotes the strength of the mean color 
field The moving particles generate a color current 
= J^t“/ 2 which forms the source term of the inho¬ 
mogeneous Yang-Mills equations 
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for the mean color field. These equations have recently 
been used in the weak-coupling limit (g 1) to simulate 
the effects of hard thermal loops [|^ on the dynamics of 
soft modes of a non-Abelian gauge field at finite temper¬ 
ature |p^ , |l^ . In this case, the colored particles describe 
the gauge field modes with thermal momenta, and the 
mean field describes the coherent motion of those gauge 
field modes which have a wave number k much smaller 
than the temperatu re T and are highly occupied. 

Here we use Eq. (|j.l|-^) to describe the interactions 
among the soft glue field components of two colliding 
heavy nuclei. Transverse modes are included by using 
a 3-dimensional spatial lattice. The short-distance lat¬ 
tice cut-off a separates the regime in transverse momen¬ 
tum where the dynamics of gluons is perturbative (large 
A:t) from that where perturbation theory fails (small fcx) 
and which is of interest to us here. The interaction with 
the mean color field allows for an exchange of an arbi¬ 
trary number of soft gluons. In combination with a par- 
ton cascade, the screening of the soft components of the 
gauge field by perturbative partons |1^J^ is taken into 
acco unt natu rally by the nonlinear nature of the coupled 
Eq. ( lO.lj-^ ). This will be discussed in a separate pub- 
lication~|l7[|. 

We represent the valence quarks of the two colliding 
nuclei as point particles moving in the space-time con¬ 
tinuum, and interacting with a classical gauge field de¬ 
fined on a spatial lattice in continuous time. Here we 
neglect the collision integrals describing hard interac¬ 
tions between the particles. In the spirit of the statis¬ 
tical nature of the transport theory, we split each quark 
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into a number Uq of test particles, each of which carries 
the fraction qg = Qojnq of the quark color charge Qq. 
For the gauge group SU(2) adopted here, each nucleon 
is represented by two quarks, initially carrying opposite 
color charge. A lattice version of the continuum equa¬ 
tions (Od-^) is constructed by expressing the field 
amplitudes as elements of the corresponding Lie algebra, 
i.e. S LSU(2). As in the Kogut-Susskind 

model of lattice gauge theory [Q we choose the temporal 
gauge Ao = 0 and define the following variables. 


Ux,i = ey::p{-igaiAi{x)) = 

^x^kl — ^x^k ^x-\-k,l ^x+k+l^ — k ^x+l,—l 
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Consequently, we have 
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for the electric and magnetic fields, respectively. The 
lattice constant in the s pat ial direction j is denoted by 
Qj. As one can see from (^) , the gauge field is expressed 
in terms of the link variables e SU(2), which represent 
the parallel transport of a field amplitude from a site x 
to a neighboring site (x + l) in the direction 1. We choose 
Ux,i and £x^i as the basic dynamic field variables and 
solve the following equations of motion 


^x,k (f-) — ^ g £x,k (f) £^x,k (^) 
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~^x-l,li^) {^x-l,kl{t) ~ Ux-l,l{t) |- 

Subsequently, we present results of a calculation which 
corresponds to a central collision of two Pb nuclei. We 
assume that 16 nucleons in a row collide on the collision 
axis (3-axis). In accordance with the spatial extension 
of a nucleon we choose a lattice extension of 1.2 fm into 
both transverse directions. The lattice spacing is taken 
aj = 0.3 fm in each direction, thus 4^ x 40 lattice points 
cover the volume of 8 nucleons in one row. The coverage 
of two complete nuclei requires much larger lattices and 
remains a challenge for more extended calculations in 
the future. The lattice is closed to a 3-torus and a dual 
lattice is superimposed on the original lattice in such a 
way that the lattice points are located in the centers of 
the cells of the dual lattice. The dual lattice cells are 
used to associate particles with lattice sites and thus to 
define the source terms in ( 0 )- 

To generate the initial configurations of the nuclei we 
randomly distribute color charged massless particles over 
the volume such that each lattice cell is occupied by an 
even number ni, of particles. The total initial charge is 
zero in each box corresponding to a neutral charge distri¬ 
bution. Momenta with opposite but random directions 
and Boltzmann distributed absolute values are assigned 
to each pair of particles. The equilibrium configuration 


of a single nucleu s in its rest f r ame is obtained through 
the evolution of (OT-^, 3^, 3J) over a long period 
of time starting from the initial fields £x,k{t = 0 ) = 0 , 
£lx,kit = 0) = I 2 . As shown in Table 1, we find ratios 
between field energy Wf and total energy ITtot which de¬ 
scend smoothly from Wf/Wtot = 0.45 to IFf/ITtot = 0.14 
for particle densities rib = 2 to Ub = 10. In Table 1 we 
also demonstrate that the ratio does not depend on a, g, 
and Qo- The ratio is also independent of the time step (if 
At < a/10). For the large final times ~ 10^, the value 
of Wf/Wtot has saturated, except for very small cou¬ 
plings, such as g = 0.3, where the equilibrium is reached 
only after longer times. 

The simulation of a collision requires a Lorentz boost 
of each nucleus into the center of velocity frame of both 
nuclei. In the example presented below, the kinetic en¬ 
ergy is 100 GeV/u for which 7 = 106.5. Both nuclei are 
mapped into their initial positions on a large Lorentz- 
contracted lattice right after the boost of particle co¬ 
ordinates Xi,pi and field amplitudes. We have investi¬ 
gated two ways to perform the boost. In one method we 
thermalise the nuclei first and apply the Lorentz trans¬ 
formation afterwards. This method is complicated and 
requires a continuation of the time-evolution throughout 
the boost. Furthermore, the lattice dispersion relation is 
not Lorentz invariant and therefore it does not lead to 
a stable translation of the nuclei on the lattice after the 
gjooost. 

More successful is a boost of the particle coordinates 
at time t = 0 for both nuclei, wherafter the fields 
cire generated from the initial conditions £x,k{t = 0 ) = 0 
and Ux,k{t = 0) = I 2 . This method leads to quasi-stable 
configurations propagating over distances which extend 
over several thousand lattice points in the longitudinal 
direction. However, it requires large initial distances be¬ 
tween the nuclei because an equilibrium between particles 
and fields has to be established before the collision. For 
nucleons with the correct mass (m„ = 939 MeV) about 
50% of the total energy are carried by glue fields. For the 
above chosen lattice constants each nucleon is covered by 
80 lattice cells and thus, according to the results in Ta¬ 
ble 1 , each quark must be split into riq = 160 particles to 
obtain a ratio close to Wf/Wtot = 0.5 at equilibrium. In 
order to reproduce the correct nucleon mass the particle 
momenta have to be smaller by a factor t/riq on aver¬ 
age in comparison to the momenta of the quarks. The 
field energy associated with a single link can then occa¬ 
sionally exceed the kinetic energy of a particle resulting 
in the emission of particles due to backwards scattering 
from the link. This corresponds to the (unphysical) emis¬ 
sion of nucleons from the moving nuclei, which we want 
to avoid. 

As long as we are only interested in the dynamics of 
the glue fields we may increase the momenta of the par¬ 
ticles such that the particle distributions of the nuclei 
remain stable throughout the whole collision. Therefore, 
in order to prepare a stable initial state for the collision, 
we start (before the Lorentz boost) with an initial condi- 
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tion in which the momenta of the particles are generated 
at a fictitious temperature of Tp = 100 GeV. The di¬ 
rections of the momenta are randomly distributed into 
the forward hemisphere such that pairs of particles have 
opposite transverse momenta. We found that the final re¬ 
sults are quite independent of the initial condition. These 
initial conditions, however, contain much higher total en¬ 
ergy than Wtot and are unphysical. The fact that 
the temperature of the gluon field T ^ Tp is limited by 
Tmax = 1-3 GeV due to the lattice cutoff prevents the 
field energies from becoming too large. Furthermore, we 
increase the value for nt to keep the ratio Wf/Wtot small 
while IFtot ^ bFtot • For nt = 6 choosen, an initial 
separation do of more than 12000 lattice spacings is nec¬ 
essary to approach a ratio of about Wf/Wtot = 0.5. do 
is Lorentz contracted to 33.8 fm in the center of veloc¬ 
ity frame. According to the above fixed extensions into 
transverse directions and size of the nuclei we use a lat¬ 
tice that comprises a total of 4^ x 12100 points. While the 
nuclei translate on the lattice, we adiabatically increase 
the coupling constant with a rate of Ag/At < (5000a)“^. 
After g has reached its final value {g = 2.0), we continue 
the time evolution until Wf/Wtot ~ 0.5. 



FIG. 1. Transverse and longitudinal T-field energy densi¬ 
ties at various times before the collision. 


In the following, we call the direction j = 3 which is par¬ 
allel to the collision axis the “longitudinal” direction and 
T-field amplitudes associated with longitudinal links are 
called the “longitudinal fields” Ey. The fields Et in trans¬ 
verse directions are called “transverse”. Figures la,b 
display the field energy densities w[^\z) = Tr{£i£i) + 
Tr(£' 2 ^ 2 ), (see Fig. la) and w[^\z) = Tr{£o£o) (see Fig. 
lb) for both, transverse and longitudinal compontents of 
the T-field every 10000 time steps. A time-step width 


of 0.03 fm has been used. We verified, that the corre¬ 
sponding i?-field energy densities (z) are identical 
with w) (z). After 1000 time steps (not shown in the 
figure) the transverse T-field energy w/'®^ is by about 
a factor 7/2 larger than the longitudinal T-field energy 
W{ '. This ratio increases to 7 at 50000 time steps and 
remains constant afterwards. 

As Fig. la shows, two configurations representing two 
nuclei propagate remarkably stably (nucleus 1 from left 
to right and nucleus 2 from right to left) over long dis¬ 
tances on the lattice. In the upper left corner of Fig. lb 
at time step 40000, the particle density (P), which is de¬ 
fined as the number of particles per bin of width a /10 in 
longitudinal direction, is superimposed on the transverse 
(T) and longitudinal (L) field energy densities for com¬ 
parison. The particle density defines the extension of the 
nucleus in the longitudinal direction. The mismatch be¬ 
tween particle density and field energy density decreases 
when the density of lattice points covering the nucleus 
in the longitudinal direction is increased. This requires a 
shift of the cutoff (here at kc = 2.0 GeV for a = 0.3 fm) 
to higher momenta. 

The tails of the field energy densities result not only 
from lattice dispersion effects. The lattice is Lorentz con¬ 
tracted in the longitudinal direction and particles have 
large longitudinal momentum components. Therefore, 
particles primarily transfer energy into longitudinal links. 
The nonlinearity of the Yang Mills equations provides 
a mechanism transferring energy from longitudinal into 
transverse degrees of freedom. Field amplitudes which 
are left behind on the longitudinal links are cancelled by 
following particles of rotated color charge. Fig. lb shows 
that this cancellation does not work perfectly, and small, 
but long tails remain behind the nuclei. The amplitude 
of these tails decreases for increasing rib and do ■ Another 
advantage of choosing rib large is to avoid artificial exci¬ 
tation of high frequency modes caused by an interplay of 
lattice discretization and point-like charges. 

Fig. 2a and Fig. 2b display w[^\z) and wl^\z) at 
4000, 8000, 12000 and 16000 time steps after the colli¬ 
sion. The particle densities are superimposed to indicate 
the position and original extension of the nuclei. As com¬ 
pared to Fig. lb (upper left window) the width of the 
distributions wl (z) moving with the particles is con¬ 
siderably increased. At the begining of the collision we 
observe a kink in Wf(t), which increases remarkably in 
the collision. As we see from Fig. 2a, a considerable frac¬ 
tion of the transverse field energy is deposited around the 
center of collision between the receding nuclei. At time 
step 76000, about 30.1 GeV of the total transverse elec¬ 
tric field energy iv/= 490.7 GeV (at time step 60000) 
are left in the spatial region between 13.9 fm and 19.9 fm. 
Consequently, the energy density in the overlap region 
with the particles is reduced allowing an enhanced trans¬ 
fer of energy from the particles into the fields. This re¬ 
sults in the kink in W{{t). 
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Fig. 2b displays the corresponding longitudinal energy 
densities. The calculation shows that within short times 
(At = 0.5 fm) after the collision about 26% (7.25 GeV) 
of the total field energy (28.03 GeV) left in the region 
between 16.4 fm and 17.4 fm are carried by longitudinal 

/ Tp\ / 

components of the FI-held. The ratio wf /w^ Mn this 
region is 0.7 which is large compared to the above men¬ 
tioned factor 1/7. Since the color charge and color cur¬ 
rent density is zero, the linearized equations (0.4) are 
homogeneous in this region and allow only for solutions 
with non-zero amplitudes into transverse directions in 
relation to the direction of their energy how. Figure 2b 
displays the fraction with momenta pointing into trans¬ 
verse direction relative to th e co llison axis. A calculation 
with the right-hand side of (0.2) set to zero yields prac¬ 
tically the same results, which hnds its explanation in a 
transfer of energy from the propagating transverse helds 
to the longitudinal hel ds due to the non-linear coupling 
between Et and Ei in (^^). When the two nuclei overlap 
during the collision, the helds are superposed. As a re¬ 
sult of the nonlinear terms in (^^), which act as source 
terms for the longitudinal helds, the amplitudes E\ start 
to grow rapidly. This mechanism describes the scattering 
of soft gluons. 



longitudinal coordinate [fm] 

FIG. 2. Transverse and longitudinal A-field energy densi¬ 
ties after the collision. 


The distribution w) (z) increases sharply for collision 
times larger 1.0 fm. When the two nuclei pass through 
each other, the particles experience the combined helds 
of both nuclei. The corresponding changes of the held 
amplitudes enter into the r.h.s. of (0.3) and modify 
the orientation of the color charge vectors ^{t) during 
the rather short overlap time of 0.11 fm. This induces 
net color charge currents resulting in radiation of glu¬ 


ons. Since the longitudinal momenta of the particles are 
large compared with their transverse components, color 
charges move essentially parallel to longitudinal links and 
induce electric helds on these links. Before the collision, 
the charge of the following particles was polarized such 
that these helds were cancelled. After the collison this is 
no longer true. A continuation of the time evolution has 
shown that the longituinal energy density of the helds 
remains as large as in Fig. 2b and decreases slowly af¬ 
ter time step 100000 (28.17fm). These longitudinal helds 
possess only transverse momenta and result in gluon ra¬ 
diation into transverse directions. 

The energy distribution in Fig. 2a and in particular 
the hump at the ce nter of collision result from non-linear 
properties of Eq. (^) and exhibit the gluon-gluon in¬ 
teraction. We repeated the calculation for smaller initial 
distances do- It turned out that there is a threshold be¬ 
tween 7000a and 9000a below of which the hump doesn’t 

/ ipN / r'\ 

appear and ic) (z) < w) {z) between the receding nu¬ 
clei. Further, (z) is almost constant as a function of 
z in this case. As can be seen from Fig. 2a and 2b, above 
this threshold, we hnd w) (z) > wl (z) in a small re¬ 
gion around the center of collision. The slow decay of the 
hump indicates that essentially low frequent modes are 
excited in the collision of the fields. The transverse fields 
have momenta into longitudinal directions and thus de¬ 
scribe longitudinal gluon radiation. A comparison of the 
transverse field energy deposit with the longitudinal field 
energy deposit for large times {t/a = 76000) shows that 
the radiation of soft gluons into longitudinal directions 
amounts to about 30% of the intensity into transverse 
directions. 
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a [fm] 

9 

<3o 

rib 

in 

[a/10] 

Wt/Wtot 

1.0 

1.0 

2.0 

2 

1 


0.45 

1.0 

1.0 

2.0 

4 

1 

• 10“^ 

0.30 

1.0 

1.0 

2.0 

6 

1 

• 10“* 

0.26 

1.0 

1.0 

2.0 

10 

1 

• 10“* 

0.14 

0.1 

1.0 

2.0 

2 

1 


0.45 

0.1 

1.0 

2.0 

4 

1 

• lO'^ 

0.30 

0.3 

1.0 

2.0 

2 

1 


0.45 

0.3 

2.0 

2.0 

2 

1 

• 10“* 

0.46 

0.3 

3.0 

2.0 

2 

1 

• 10“^ 

0.46 

0.3 

0.5 

2.0 

2 

1 


0.37 

0.3 

0.5 

2.0 

2 

2 

• 10“* 

0.44 

0.3 

0.3 

2.0 

2 

1 


0.19 

0.3 

0.3 

2.0 

2 

3 

• 10“* 

0.36 

0.3 

0.3 

2.0 

2 

6 

• 10“^ 

0.45 

0.3 

2.0 

1.0 

2 

1 


0.46 

0.3 

2.0 

0.5 

2 

1 

• 10“* 

0.46 

0.3 

2.0 

0.3 

2 

1 

■ 10“* 

0.46 


TABLE 1. The ratio WtjWtot is listed in dependence on 
the parameters a, g, Qo, rih for large times tn when the con¬ 
figuration has evolved close to equilibrium. 


















